Nonlinear Transport through ultra-narrow Zigzag Graphene Nanoribbons: 
non-equilibrium charge and bond currents 



Hosein Cheraghchi 

School of Physics, Damghan University, P. O. Box: 36715- 364, Damghan, IRAN^ 

(Dated: November 14, 2012) 

The electronic nonlinear transport through ultra narrow graphene nanoribbons (sub-lOnm) is 
studied. A stable region of negative differential resistance (NDR) appears in the I-V characteristic 
curve of odd zigzag graphene nanoribbons (ZGNRs) in both positive and negative polarity. This 
NDR originates from a transport gap inducing by a selection rule which blocks electron transition 
between disconnecting energy bands of ZGNR. Based on this transition rule, on/off ratio of the 
current increases exponentially with the ribbon length up to 10 5 . In addition, charging effects and 
also spatial distribution of bond currents was studied by using non-equilibrium Green's function 
formalism in the presence of e-e interaction at a mean field level. On the other hand, we also 
performed an ab initio density functional theory calculation of transmission through a passivated 
graphene nanoribbon to demonstrate robustness of the transport gap against hydrogen termination 
of the zigzag edges. 
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I. INTRODUCTION 



Graphene is a two-dimensional carbon which has 
been recently explored^. Experiments in graphene-based 
devices^ have shown the possibility of controlling their 
electrical properties by the edge structure and also ap- 
plication of an external gate voltage. Nowadays, by us- 
ing a chemical methocP, it is possible to fabricate field 
effect transistor graphene nanoribbons (GNRs) with ul- 
tra narrow widths and smooth edges which are possi- 
bly well defined zigzag or armchair edge structures. A 
useful transport gap is opened in sub-lOnm GNRs at 
room temperature which results in a high on/off current 
switching up to 10 6 . The origin of the transport gap can 
be understood by two factors: suppression of transport 
due to edge disorder leading to Anderson localization, 
confinement along transverse directiorP^. Additional to 
these factors, in ultra narrow zigzag graphene nanorib- 
bons (ZGNR), flow of current is also blocked by the sym- 
metric selection rules^E^. 

Negative differential resistance^ in nanoelectronic de- 
vices has been also observed in the metallic nanotubc 
junctions^ and in the case of potential barriers in 2D 
graphene sheets^. 

In this paper, to shed light on the experimental work 
in Ref.|3]j we study non- linear transport in odd zigzag 
graphene nanoribbons by using non-equilibrium Green's 
function formalism (NEGF). It is shown that a stable 
NDR against electrostatic interaction up to 10 5 is ap- 
peared in ultra narrow ZGNRs around ±li in both posi- 
tive and negative polarity, where t is hopping integral be- 
tween carbon atoms. On-off ratio of this NDR increases 
exponentially with the ribbon length. The NDR phe- 
nomenon occurs in far-from-equilibrium regime. 

Although parity is not conserved in odd ZGNRs, cur- 
rent reduction is induced by an other selection rule in 
which transition between disconnecting band groups is 
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FIG. 1: Odd zigzag graphene nanoribbon considered as a cen- 
tral region attached to two electrodes. Transverse zigzag car- 
bon chains and unit cells inside the central region are labeled 
by N and M. 



forbidden (Cresti et o/P). The edges of graphene ribbons 
can also simply absorb some chemical compounds^. It 
is demonstrated that this NDR is not much sensitive to 
the asymmetry of ribbons, while transport gap in even 
ZGNRs which is based on the reflection symmetry, is 
completely destroyed due to the edge disorder. The above 
mentioned selection rule is valid for both even and odd 
ZGNRs, while parity selective tunneling just belongs to 
even ZGNRs. This transport gap is robust against hy- 
drogen termination. Additional to the model calculation, 
by using TranSIESTA codtP, we performed an ab ini- 
tio density functional theory calculation of transmission 
through graphene nanoribbon with hydrogen termination 
of the zigzag edges. Results of first principle calculation 
confirms the transport gap originating from transition 
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rule between disconnected bands. 

This NDR is stable against the presence of the electron- 
electron (e-e) interaction at a mean-field approximation. 
It is concluded from self-consistent charge and potential 
profiles so that external potential is screened by charge 
redistribution around the contact junctions. So potential 
profile deep inside the ribbon remains flat. Furthermore, 
the e-e interaction increases on-off ratio of the current. 
Moreover, in low and high applied biases, we study spa- 
tial profile of local currents in the presence of e-e interac- 
tion which contrasts with non-equilibrium charge profile. 

The paper is organized as the following sections. Af- 
ter an introduction, we explain Hamiltonian and a short 
review of NEGF formalism in section II. In section III, 
we present our results containing NDR in I-V curve and 
charging effects in high applied bias and spatial bond 
currents. Before to end the paper, we present the results 
of ab initio calculated by using Transiesta code. Finally, 
we conclude our results. 



II. FORMALISM 

The single electron Hamiltonian of the molecule is de- 
fined as: 




FIG. 2: Current-voltage characteristic curve of an odd zigzag 
graphene nanoribbon with N = 5 (zigzag chains). I-V curves 
are compared for two ribbon lengths; M = 5 and 10. For the 
case of (M, N) = (5, 5) and for comparison purpose, I-V curve 
is also plotted in the presence of electron-electron interaction 
(U) and asymmetric factor(e a ). 
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where c\ and Ci are the electron creation and annihila- 
tion operators, respectively. 27V x M and t are number of 
atomic sites and hopping energy between nearest neigh- 
bor atoms. One n orbital is considered per each site for 
graphene as a planar 2D system. Without losing any gen- 
erality, we set onsite energies (e^) of all sites to be equal to 
zero. All energies and voltages are in unit t = 2.7eV. The 
Fermi energy of undoped graphene nanoribbon is at half- 
filling Ep = due to perfect electron-hole symmetry. By 
applying a source-drain bias, the site energies are shifted 
by a linear potential distribution along the molecule itf 
which is the solution of the Laplace equation. The ap- 
plied bias V is divided symmetrically on the left and 
right electrodes as — V/2 and V/2, respectively. Vtj is 
the electrostatic Green's function and Srii = i%i — n® is 
the change in the self-consistent charge rii from its value 
in zero source-drain voltage n®. This third term is the 
direct Coulomb interaction created by the source-drain 
bias-induced charges. It includes the Hartree term which 
is the solution of Poisson equation and locates on the 
diagonal terms of Hamiltonian. For calculating the elec- 
trostatic Green's function, we have used the method ex- 
plained in Ref. |17j and its appendix. For self-containing, 
we present a very brief review of the NEGF formalism. 
Charge density in non-equilibrium situation is calculated 
by [— iG < ] as the occupation number in the presence of 
the two electrodes with an applied source-drain biaJ^. 



T C Ef ~t 



Im[G r {E)] u dE + nl 



(2) 



where non-equilibrium part of charge can be calculated 
by the following integral, 
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[-iG<(E)} u dE (3) 



where within a one-particle theory 



iG< = G r (T L f L +T R f R )G a 



(4) 



Here Jl/r is the Fermi-Dirac distribution function of 
electrodes and G r /a is the retarded/advanced Green's 
function defining as the following: 



G r/a = p ± _ H { n y _ ^r/a _ £ « 



(5) 



and r is the escaping rate of electrons to the electrodes 
which is related to the self-energies as T p — — E£] 
with p = £/iJlSH20|. Here 77 — ^ + . Solving equa- 
tions [2] and [5] self-consistently results in a self consistent 
charge and Green's functions. Finally, the current pass- 
ing through the molecule is calculated by the Landauer 
formula for zero temperature^! which is valid for coherent 
transport. 
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FIG. 3: Energy spectrum of the left and right electrodes 
and transmission through zigzag graphene nanoribbon with 
(M,N) = (10,5) at voltage {Vsd = 1.4t > V T ). The band 
structure is divided into three groups which are called by up- 
per , centeral and lower band groups. These groups are clas- 
sified based on the bands which are connected in terms of k x . 
The bold hollow arrows show the current integration window 
which based on Eq.[6] is proportional to Vsd- The Fermi level 
is set to be as Ef = 0. The half- width of the central bands 
at k x = is called as 8 which is equal to the threshold NDR 
voltage. A is energy separation of the upper bands from the 
central bands at the Dirac points. Transport gaps AB and 
CD are equal to A — S + Vsd ■ 



9 f E F +V/2 
I(V) = -r dET(E, V) (6) 

n JBp-V/2 

where T(E, V) are the bias dependent transmission co- 
efficient. 

T = Tr[G r T R G a T L ] (7) 

Fig. ([I]) shows zigzag graphene nanoribbon with an odd 
number of zigzag chains in width (N). Central interact- 
ing region, left and right electrodes divide the ribbon into 
three regions. Odd ZGNRs has a bisection plane which 
its intersection with the ribbon has been shown as the 
x-axis in Fig. ([I]) . Although this plane bisects the ribbon, 
it is not a reflection plane. If upper half of the ribbon is 
displaced by a/2 in respect to the lower one, mirror sym- 
metry is achieved against x-axis. a is the bond length of 
C-C. In other words, part of the wave function which is 
a functional of y variable, is an eigenvector of the parity 
operator^! 

III. RESULTS 

Current-voltage characteristic curve of an odd ZGNR 
with 5 zigzag chains (N — 5) is shown in Fig.|2]). Lower 



than the external bias 1.2t, current increases linearly with 
the applied bias as an Ohmic device. After a threshold 
voltage (1.2t), NDR occurs at both positive and negative 
polarity. 

The origin of NDR seen in odd ZGNRs is inter- 
preted by analyzing their energy spectrum accompanied 
to transmission curve. Fig. ^ shows energy spectrum 
E(n, k x ) and transmission through 5-ZGNR at V = lAt. 
n counts bands from the bottom (n = 1) to the top of 
band structure (n = 27V) and k x is the longitudinal mo- 
mentum. 

In odd ZGNRs, parity has noncommutative relation 
with the Hamiltonian. Therefore, parity has no con- 
servation and consequently transmission is not blocked 
by the parity selection rule, while parity conservation in 
even zigzag nanoribbons opens transmission gap around 
Fermi levef^. In the range of BC of Fig.([3]), there exists 
one conducting channel which results in the unity trans- 
mission around the Fermi level. In this range, electrons 
which occupy — k states of the lower band of the cen- 
tral band group (the dashed blue bands) belonging to the 
right electrode are injected into the unoccupied — k states 
of the upper band of the central group belonging to the 
left electrode. So at low biases, current increases linearly 
with the bias. This single-channel transport around the 
Fermi level remains unchanged even for high voltages. 
However, for voltages greater than the NDR threshold 
voltage (V > Vt), blocked regions marked by the ranges 
of AB and CD comes into the current integration win- 
dow. The current integration window, based on Eq. [6] 
and Ef = 0, is proportional to Vsd and is shown with 
the bold hollow arrows in Fig. ([3]) . Therefore, when the 
source-drain applied bias increases, current begins to de- 
crease. 

Blocked regions (AB and CD) arise from a selection 
rule which increases back scattering in the lengthy rib- 
bons. According to this rule, electron transition between 
those bands which are disconnected from the view point 
of longitudinal momentum, decreases exponentially with 
the length. Topology of zigzag graphene ribbons di- 
vides the band structure into the three different groups^. 
Three groups which are called by upper (n=l,2,3,4), cen- 
teral (n=5,6) and lower (n=7,8,9,10) groups, are classi- 
fied based on the bands which are connected in terms 
of k x . Variation of the momentum of electrons passing 
through the system strongly depends on smoothness or 
sharpness of the potential. The transition probability of 
electrons with the k state injected by the right electrode 
into the q state as an empty state in the left electrode is 
proportional to the Fourier transform of the longitudinal 
voltage^ as < ip(k) \ V | i/j(q) >oc V(k — q). So smooth 
variation of the potential in longer ribbons results in a 
small momentum variation of electron. Consequently, a 
smooth potential in the longitudinal direction can just 
scatter electron into those states belonging to the energy 
bands which are connected from the point of momentum 

viCT^SI. 

Bias dependent of transmission leads us to plot a con- 
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FIG. 4: Contour plot of transmission versus energy and ap- 
plied bias for graphene nanoribbon with 5 zigzag chains in 
width and 10 unit cells in length. Dark solid lines show 
the current integration window and white dashed lines show 
the forbidden region for electronic transition from the band 
groups of upper or lower bands to the bands of central group. 



tour of transmission in plane of energy and applied bias 
as shown in Fig.Q. In this figure, the solid lines show the 
current integration window. Moreover, the dashed lines 
clarify the regions AB and CD corresponding to trans- 
port gaps shown in Fig.([3|. The intersection of these 
blocked regions and also the current integration window 
would be around 1.2t which is the threshold voltage Vr 
for the current reduction. 

Band structure analyzing demonstrates that the 
threshold voltage is equal to the half-width of the cen- 
tral bands at k x = as Vr = 6 = [E(N + l,k x = 
0) - E(N, k x = 0)]/2. As shown in Fig.Q, A is energy 
separation of the upper bands from the central bands at 
the Dirac points. There is a Log-Normal behavior of S 
versus number of zigzag chains (N) such that as N — » oo, 
the NDR threshold voltage asymptotically approaches to 
the value of 0.9738 ± 0.0002t. So the NDR threshold 
voltage slightly decreases with the ribbon width. 

Analyzing transport gaps appeared in the band struc- 
ture shows that they are equal to A — 6 + Vsd where 
A oc N . Since 8 approaches to a constant values when 
N — > oo, in a given voltage, transport gap is disappeared 
for N > 30 which is nearly equivalent to lOnm. 

The other effect which enhances performance of this 
electronic switch, is the ribbon length. Fig. ^ shows 
an increase of on/off ratio with the ribbon length. More- 
over, NDR region (V f / — V on ) occurs in a more extended 
range of the I-V curve. Exponential decay of transmis- 
sion with the length in the gap regions develops quality of 
switching. It is shown in Fig.flBl that I on /I ff increases 
exponentially with the ribbon length. 

Another parameter which affects I-V curve, is the e- 
e interaction. In this case, interaction intensifies NDR 
effects such that off current reduces in compared with 
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FIG. 5: Exponential increase of on/off current ratio as a func- 
tion of ribbon length (M). 



non-interacting system and also the region of NDR be- 
comes more extended. 

If one of the ribbon edges is doped by small impurity 
such as e a = O.li, because of a band gap which is in- 
duced by edge impurity at the band center, current at 
low biases decreases. Fig. ^ represents that even with 
the presence of edge impurity, still the region containing 
NDR exists. However, asymmetry decreases on/off ratio 
of the current. Furthermore, asymmetric ZGNR behaves 
as a semiconductor while symmetric ZGNRs behave as 
an Ohmic devices^. The effect of asymmetry on NDR 
competes with the ribbon length. Since asymmetry can 
not be ignored in experiment, longer ribbons are in favor 
of keeping NDR in the I-V curve. 

On the other hand, additional to edge impurity, this 
asymmetry can be assigned to a sublattice symmetry 
breaking induced by spontaneous ferromagnetic spin or- 
dering of the electrons localized at the zigzag edges^. 
In fact, the border atoms at the two opposite zigzag 
edges belong to different sublattices. So spin orienta- 
tion along the edges induces different magnetic poten- 
tials at the edges. As a result, a small band gap is opened 
around Fermi level which depending on the ribbon width, 
is about 0.15 eV. The asymmetry which we have con- 
sidered is about 0.3 eV which is stronger than the gap 
opened by spin-polarization of the edges. We can con- 
clude that spin-polarization along the zigzag edges can 
not affect emerging of this NDR phenomenon^. 

Charging effects: To understand why details of elec- 
trostatic potential do not affect the emergence of NDR, 
we compare transmission in the presence and absence of 
e-e interaction. Figure ^ compares transmission curves 
for voltages 1.2t at the NDR threshold Vr and a larger 
value such as 1.5i. It is clear that for voltages lower 
than the NDR threshold, transmission in the conducting 
region BC is robust against e-e interaction. However, 
for voltages V > Vr, transmission of conducting region 
BC reduces with e-e interaction. This is the reason for 
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Ribbon Axis(A) 







~- > 









FIG. 6: The effect of efectron-electron interaction in the 
mean-fieid ievei on transmission curves for two voltages a) 
NDR threshold voltage which is equal to 1.2t, b) a voltage 
above the NDR threshold such as 1.5t. 



reducing off-current in the presence of e-e interaction. 
Furthermore, transmission curve in the blocked regions 
CD and AB exhibits a little enhancement with the e-e 
interaction for the whole range of voltages. 

To elucidate physics behind to the robustness of 
NDR phenomenon against e-e interaction, we investi- 
gate self-consistent potential profiles which are shown in 
Fig.(7ja,b). In fact, external potential is well screened 
by charge redistribution so that potential mostly drops 
at the contact junctions. This fact is represented in 
Fig. ([7]). So, electrostatic potential of atoms located far 
from contact junctions (deep inside the ribbon) remains 
nearly flat. In more details, depletion of charge close to 
the source electrode reduces the source potential as USn. 
Therefore, source potential can not penetrate inside the 
central portion of the system. On the opposite side of 
system, charge is accumulated close to the drain elec- 
trode such that dropping of the potential is weakened. 
It is useful to investigate bias dependence of total trans- 
ferred charge from/into the system. Fig.([7]c) shows that 
transferred charge in compared with its neutrality point 
value (no with zero source-drain applied bias) decreases 
for voltages higher than the NDR threshold, while it re- 
mains unchanged for voltages V < Vr- In other words, 
for V < Vr, in- flowing and out-flowing charges are bal- 
anced with each other, while for V > Vr, system becomes 
empty of charge. It can be seen that the threshold voltage 
for depletion of charge is correlated with the NDR thresh- 
old voltage. This behavior is closely similar to what is 
seen in graphene nano-junctions 17 . As a consequence, 
at high voltages, due to accumulation of charge close to 
the drain electrode, external potential is always screened 
in the middle of the graphene strip. Screening effect in- 
duces sharp variation of electrostatic potential at the con- 
tacts which enhances transition probability between dis- 
connected band groups. Therefore, transmission in the 



FIG. 7: a) Self-consistent charge profile at 1.5t. Charge com- 
pares with its neutrality point (no) when there is no applied 
bias. Transferred charge is equal to (n — no). Electrostatic po- 
tential profile is calculated according this transferred charge, 
b) averaged electrostatic potential variation on unit cell as 
a function of the ribbon axis, c) Total transferred charge 
from/into the ribbon as a function of source-drain voltage. 



blocked regions (CD and AB) increases in compared to 
non-interacting system. However, blocked regions have 
not any contribution in the current calculation. 

Spatial distribution of bond currents: Charge conserva- 
tion based on continuity equation yields non-equilibrium 
bond charge current)^ for a non- interacting tight-binding 
Hamiltonian. 



Jij — 
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[G<AE)~G<(E)]dE 



(8) 



where sites i and j are nearest-neighbor atoms which 
their hopping integral is nonzero. Now, charge continuity 
equation is derived by using Heisenberg equation. 



e dt U * 



Jn] = 



(9) 



;th 



where j are the nearest-neighbor atomic site around 
atomic site. It can be simply proved that the integrands 
in Equations [3| and |8| are real. In fact, transpose of matrix 
—iG < in Eq^H) is equal to its conjugate. So, diagonal 
terms of (— iG>) and also the terms of (G^ — GfA are 
real. 

Bond charge current formulated in Eq(8] has been 
derived for non-interacting tight-binding Hamiltonian. 
However, by using Heisenberg equation, it can be sim- 
ply proved that this formula is still applicable for the 
Hamiltonian described in EqJT] in which Hartree inter- 
action appears on diagonal elements of the Hamiltonian. 
On the other hand, charging effects originating from elec- 
trostatic interaction (Hartree term) adjust local currents 
by means of lesser Green's function. 
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FIG. 8: Top: Spatial distribution of bonds current in three 
voltages: red arrows) 1.0 t< Vt, black arrows) 1.2t = Vt 
and blue arrows) 1.5 t> Vt- The length of arrows is propor- 
tional to the local current. The arrows located on the edges 
are shorter than the arrows located in the bulk of nanorib- 
bon. By increasing the source-drain voltage, for voltages 
lower than Vt, current increases and consequently arrows are 
longer, while higher than Vt, current decreases and so arrows 
are shorter. 

Bottom: Non-equilibrium charge profile at the threshold NDR 
voltage. Non-equilibrium charge reaches in its maximum 
value at the ribbon's edges. Moreover, in agreement with 
the continuity equation, atomic sites with larger local current 
corresponds to smaller non-equilibrium charge. 



FIG. 9: Top: Zigzag graphene nanoribbon terminated by Hy- 
drogen atoms. Bottom: First principle calculation of trans- 
mission through the above system as a function of energy for 
different values of external biases. 



that atomic sites with larger local current corresponds to 
smaller non-equilibrium charge. On the other hand, non- 
equilibrium charge increases for atomic sites with higher 
onsite energy due to applied bias (close to the source elec- 
trode), while its value tends to zero for sites close to the 
drain electrode. 



Fig.([8]top) shows spatial profile of local current at each 
site. In the center of the ribbon, magnitude of local cur- 
rent densities are larger than its value in the zigzag edges 
of the ribbon. This feature is valid for both ranges of low 
and high voltages not more than 1.5 t. The reason can be 
followed by single-channel transport through ZGNR even 
for high voltages^. One channel transport is an indirect 
result of transition rule in which transition between dis- 
connected bands is forbidden. For voltages higher than 
1.5 t, alignment of local currents along the ribbon axis is 
gradually disturbed. 

In contrast with the local current distribution, as 
shown in Fig. ([8] bottom), non-equilibrium charge reaches 
in its maximum value at the ribbon's edges. This re- 
sult is compatible with the continuity equation (Eq|9]) so 



Comparison of odd and even ZGNRs: There are some 
interesting differences between results arising from odd 
ZGNRs with those results belonging to even ZGNRs^ 
which is kind of odd-even effect. Here, NDR appears 
in voltages upper than It while in even ZGNR, NDR 
occurs for voltages lower than It. On/off ratio of the 
current in gated even ZGNRs increases as a power law 
with the function of the ribbon length while here, on/off 
ratio increases exponentially. Screening of the external 
bias by electrons of system in even ZGNRs is so stronger 
than screening effects in odd ZGNRs. As a consequence, 
the effect of electrostatic interaction on increase of the 
on/off ratio in even ZGNRs is much effective than in odd 
ZGNRs. Furthermore, transferred charge from/into the 
central portion of graphene nanoribbon depends on odd 
or even zigzag chains in width. 
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IV. HYDROGEN TERMINATED 
NANORIBBON 

To show how passivation of zigzag edges of graphcne 
nanoribbon affects transport properties, we compare the 
results of the presented model with transport proper- 
ties calculated by using TranSIESTA code^. This code 
is based on density functional theory (DFT) approach. 
We have used the following options to calculate SIESTA 
code: the generalized gradient approximation (GGA) 
with the Perdew-Burke-Ernzerhof exchange-correlation 
functional (PBEjPH, double-C plus polarization orbital 
(DZP) bases for all atoms, and Troullier-Martins norm 
conserving pseudopotentials to represent the cores, 200 
Ry real space mesh cutoff for charge density and a super- 
cell within 20 A of vacuum between the periodic graphcne 
nanoribbons. 

Fig.(|9j Top) shows a schematic view of odd zigzag 
graphene nanoribbon which is saturated by hydrogen 
atoms. Transmission around the Fermi energy is rep- 
resented in Fig. ([9j bottom). The energy scale is shifted 
so that the Fermi energy of the system, if there is no bias 
voltage, is zero (Ep = 0). When there is a finite bias, the 
Fermi energy of the left electrode is placed at V/2 and 
the right electrode at —V/2. The occupation number on 
atoms in the center part of the system is determined by 

As it is clear of Fig.([9j bottom), there exists one trans- 
mitting channel for each spin around Fermi energy level. 
This result is in complete agreement with the band struc- 
ture analysis shown in section. III. There exists only one 
transmitting channel in the region BC presented in fig- 
ures [3] and |4j The transport gap in Fig. ([9j bottom) 
begins to open at voltage around 1.2eV ~ 0.5t and for 
energy about 2.7eV ~ t. These points correspond to the 
points marked by cross sign in Fig|4]in which the trans- 
port gap begins to open as the applied bias increases. As 



a result, we demonstrate that the transport gap which 
is responsible for emerging the NDR phenomena could 
also be opened in graphene nanoribbons with passivated 
zigzag edges. 

V. CONCLUSION 

As a conclusion, based on a model calculation of non- 
equilibrium Green's function formalism, we found that 
there exists a region of negative differential resistance 
in I-V curve of ultra narrow (lower than 10nm) zigzag 
graphene nanoribbons with odd number of zigzag chains 
in width. This NDR is induced by a transport gap which 
originates from electronic transition between disconnect- 
ing bands of energy from the view point of longitudinal 
momentum. On/off ratio of the current exponentially in- 
creases up to 10 5 as a function the ribbon length which 
proposes possibility of manipulation of odd ZGNRs as 
high quality switch in nanoelectronic based on graphene 
nanoribbons. In addition, e-e interaction enhances on-off 
ratio of the current which originates from a flat electro- 
static potential deep inside the ribbon due to screening of 
the external bias by electrons close to the junctions. By 
using continuity equation, spatial profile of local currents 
is calculated in the presence of Hartree electron interac- 
tions. In both high and low biases, local current reaches 
to its maximum values in the center of the ribbon while 
in contrast with the local current profile, non-equilibrium 
charge has its maximum values at the edges of the rib- 
bon. Furthermore, this NDR is not much sensitive to 
the edge asymmetry. So emerging of this NDR is robust 
against spin orientation along the edges. 

By using ab initio density functional theory, we also 
show that the transport gap which is responsible for 
emerging NDR, exits in passivated graphene nanoribbons 
by Hydrogen atoms. 
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